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Abstract 

Observations regarding rotating cavitation and 
cavitation surge experienced during the 
development of the Fastrac engine turbopumD are 
discussed. Detailed observations acouired from 
-e analysis of both water flow and liquid oxvgen 
ast data are offerea in this paper. Scaling ana 
general comparison of rotating cavitation oetween 
water rlow and liquid oxygen testing are 
discussed. Complex data features linking the 
ocanzea rotating cavitation mechanism of the 
nducer to system surge comDonents are 
jescriDea in detail. Finally a aescriDtion or a 
jmoea-oarameter nvarauiic system moaei 
,:eveioDea to oetter understand ooservea data is 
given. 

Nomenclature 

•j synchronous oscillation 

otating cavitation osculation 

C dimensional pump compliance 

K dimensionless pump compliance 

•y dimensional mass flow gain factor 

.1 timensiomess mass flow gain factor 

,ss action specific speed 
met tip flow coefficient 
t nlet cavitation number 

,t. mducer tip cavitation numoer 

■p stage head coefficient 


j, ntroduction 

During the deveiooment of the nigh speed 
:_irooDump for the 60.000 pound thrust Fastrac 
'ocket engine, complex unsteady flow ana shaft 


vibration due to rotating cavitation have been 
observed. The bipropellant pump (Figure 1 ) with 
both a single-stage LOX (shown with axial inlet) 
and RP inducer-impeller sharing the same drive 
shaft exhibits independent cavitation induced 
suDer-syncnronous oscillations concurrently during 
"ominal ooeration. 



Figure i. Fastrac TurboDumc Assembly 


"his cocument attempts to summarize maior 
observations regarding the Fastrac LOX 
turbopump rotating cavitation anomaly and relay 
them to a turbomachinery community that must 
frequently contend with hydraulic system 
nstabilities induced by cavitating pumos. 
Dbservations nave been derived from LOX duitid 
super-scale model water flow testing and 
turbopump component level hotfire test. 
Characterization and abatement of the cavitation 
instability continue to be critical issues in the 
development and certification of the Fastrac 
engine due to leading edge blade deformations 
found in both the initial aluminum water flow test 
article and subsequently in a component level test 
inducer. Figures 2 and 3 show the damaged 
leading edge regions in the water flow aluminum 
and the LOX K-Monel inducers respectively. 
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latitude offered by the water rig testing has been 
crucial to the characterization of the oscillation. 



Figure 2. Failed Aluminum Inducer 



gura 2. LOX inducer with Slight Leading tdge 

Deformation 

Unsteady 'low analysis of the Fastrac engine 
_OX turoooumc nas benefited from multipie 
sources or exoenmental/test aata of the comoiex 
nstaDiutv. Water flow rig testing of the LOX DumD. 
•’itiatea to experimentally verify suction 
oertormance. nas provided vaiuaole insight into 
the anomaly and produced the first evidence of the 
severity of the rotating cavitation with the failure of 
the original aluminum superscale LOX 
nducer-impeiler Figure 4 is a loint 
- me-freouencv ~aooinq of the oscillation taKen 
-om water how 'est or the lQX noucer at the 
'ated flow condition io = 0.135) ano a shaft soeeo 
or 5000 RPM. 

The inducer inlet fluctuating pressure data 
shows the development of the rotating cavitation 
oscillation, denoted as ‘re" in the figure, during an 
nlet pressure ramp test (suction pressure 
decrease going into oagei at constant 
•:ow-to-soeed ratio. In-depth discussion of the 
spectral trending will be offered in subsequent 
sections of this overview. The experimental 
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Figure 4. Waterfall Plot of Inducer Inlet 

Fluctuating Pressure from Water Flow Test 

Several of the features of the Fastrac 
turbopump rotating cavitation have been 
documented by many researchers involved in the 
development of other high speed rocket engine 
turbomachmery. ' 3 Tsu|imoto. et al.. provided an 
excellent overview of inducer rotating disturbances 
encountered during the development of the LE-7 
LOX turbopump. A maioritv of the spectral 
components shown in Figure 4 have counterparts 
described in the LE-7 cavitation oscillation 
overview. 

Rotating cavitation nas proven to be more 
than just a nuisance to the Fastrac engine 
development program. Along with deforming 
inducer leading edges, pump cavitation nas 
hampered characterization of the turbopump 
vibration environment due to saturation of 
accelerometers and signal conditioning. Excitation 
of both LOX ano RP feed-svstem nydrauiic modes 
py the active inducers nas aiso complicated 
engine calibration efforts by hindering accurate 
propellant flow measurement. 

II. Rotating Cavitation and Cavitation Surge Data 
Observations 


Rotating Cavitation 

As noted previously, water flow testing offered 
several experimental opportunities not available in 
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the LOX environment. This included the 
accessibility of high frequency instrumentation 
locations and the ability to easily manipulate pump 
inlet hydrodynamic condition, i.e. Nss (suction 
specific speed) and o, in the water flow loop 
(Figure 5). Flush mounted fluctuating pressures at 
the inducer inlet plane approximately 0.3 inducer 
diameters upstream of the blade leading edge 


(instrumentation plane 1000, Figure 5) were used 
to establish the rotational character of the 
oscillations shown in the joint time-frequency plot 
of Figure 4. Hannover plotting of inducer inlet 
fluctuating pressure relative phase information 
was used in determining cell number and 
propagation direction of the major spectral 
components.' 



7 gure 5. MSFC_ !naucer_ • esr Loco and Fastrac lOX Inaucer Model 


As exoected. the peaks denoted as "2N and 
N". are tnree and single ceil disturbances. 
r espectively. Cross-cnannei phasing at the 
'nducer inlet Diane confirm their lock-in un a 
rotational sense) to the three-biaoed inducer, i he 
peak denoted as rc in Figure 4. is a one-celied 
disturbance propagating in a forward direction, i.e. 
n the same direction as the rotor, at a slightly 
super-synchronous rate. No significant loss in 
pumping efficiency was noted in water flow suction 
performance mappings until the appearance or the 
synchronous oscillation, N, where stage head 
coefficient, 4 / , dropped approximately 2 percent 
from the non-cavitated (< 1 5000 Nss) water flow 
value of 0.40. 

The oeax cenotea as "fl" in Figure 4 is 
rhougnt to pe a rotating disturbance of hign ceil 
count related to backflow vortices. Phasing 
suggests that the low amplitude oscillation is an 
eight-cell disturbance with a propagation velocity 
of approximately 0.32N. The cell train propagation 
rate and high cell count compare well to rotating 
backflow disturbances characterized in 
Reference 3. However, inlet plane phase mapping 


,vas not of high enougn fidelity to confirm 
propagation direction ana the mgn ceil count of fl . 

Figure 6 shows Doth normalized rreauency 
and amplitude of 3N. rc. ana N versus Nss at the 
90% rated flow condition from both water flow and 
LOX testing. Comparison at this slightlv 
orf-nommal flow condition was selected since a 
majority of the LOX test data was acauired at the 
90% rated flow, in the figure, ciosea svmDOls 
denote normalized water data while the larger 
open symbols represent LOX. The frequency of 
the rotating cavitation component, rc, is 
normalized by synchronous frequency while all 
amplitudes are scaled by impeller tip dynamic 
pressure. Normalized freauencv for rc ranges 
! 'om 1 to t .25 and the cecreasing trend with Nss 
compares well between LOX and water test. 
Likewise, there is good comparison between 
normalized amplitudes of the rc and N 
components as well as the critical Nss at which rc 
transitions into the synchronous response. A 
normalized 3N amplitude is not shown for LOX 
since the turbooump acquired no appreciable 
steady-state operation at low Nss. i.e. less than 
12 kNss. 
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Figure 6. Normalized Frequency 3 na Amontude or Major i-astrac lGX Inducer Oscillations from Both 
~ :_OX ana Water Testing 


Cavitation Surge 

During water flow testing of the suoerscaie 
Fastrac LOX inducer model, strong pumo-facuitv 
coupling was experienced. Figure 7'sa waterfall 
plot of an upstream (approximately 1 5 inlet 
diameters from inducer) unsteady pressure. 
Amplitude scale is eguivalent to that of the inducer 
nlet waterfall plot of Figure 4 The dominant 
oscillations registered at this upstream location are 
3(rc - N) and the component denoted as 12". 
Further analysis into phasing and the nigh 
frequency content of measured fluctuating 
pressure both at the inducer inlet plane and 
elsewhere in the water flow loop revealed 
nterestmg attributes of the 3(rc - N) comDonent. 

Figure 8 is a senes or waveforms iO - ' 0 kHz 
oanowidth) across several locations over 
1 5 cycles of 3(rc - N). In sequence from top to 
bottom are upstream fluctuating pressure, inlet 
fluctuating pressure, inlet fluctuating pressure at 
an alternate circumferential location, and finally 
radial acceleration measured on the pump 
housing. Clearly, the dominance of the 3(rc - N) 
n the upstream pressure is shown. Likewise the 
complexity of the nonlinear interaction between 
the rotating cavitation component and the 
synchronous rate is very evident in the two inducer 
nlet pressures. 


While the rotation-induced delay of the higher 
‘requency rotating cavitation wavelet is easy to 
discern between the two inlet pressures, the 
underlying 3(rc - N) wave appears to be in phase. 
The inducer inlet plane transfer function of 
Figure 9 confirms this unique phasing at the 
‘requency of 3(rc — N). Relative phase at the 
oscillation freouency is zero with a coherence of 
nearly unity. Similar phasing and coherence for 
3(rc - N) were shown witn the alternate inducer 
nlet cross-channel pair. 

Before discussing possible physical 
implications of the complex component 3(rc - N), 
some observations on the transient envelope of 
the monitored fluctuating pressures are now 
offered. Figure 10 dismays a high-pass filtered 
version of Figure 8. The data were high-pass 
filtered at 2 kHz to safely exclude rc and 
synchronous activity. Of interest in this series of 
filtered waveforms is the obvious synchronization 
of wideband bursts at the rate of 3(rc - N) across 
all fluctuating pressures and monitored pump 
vibration. Envelope spectral analysis of these 
waveforms has confirmed this synchronization. 
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Figure 7. Waterfall Plot of Upstream Fluctuating 
Pressure from Water Flow Test 
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Figure 8. Multiple Location Pressure and 
Vibration Waveforms During Rotating Cavitation 
Over 15 Cycles of 3(rc - N) 
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Figure 10. Multiple Location Filtered Pressure and 
Vibration Waveforms During Rotating Cavitation 
Over 15 Cycles of 3(rc - N) 


The unique phasing is thought to reflect the 
complex physical relationship between the 
localized inducer rotating cavitation disturbance 
and the more global cavitation surge instability. A 
component similar to the 3(rc - N) has been 
documented by Tsujimoto and fellow researchers. 3 
In their effort, a cavitation surge mode was 
identified at a frequency corresponding to (re - N). 
The frequency variation of their cavitation surge at 
times lost correlation with those of the rc and N 
components. In Fastrac LOX inducer water flow 
testing, correlation of the cavitation surge mode at 
3(rc - N) with the rotating cavitation ana 
synchronous components was maintained. 


3(rc - N) 


rc 




3(rc-N) 


Figure 9. Inducer Inlet Pressure 
Cross-Channel Coherence, Phase, and Gam 
During Rotating Cavitation 


Figure 1 1. Instantaneous Frequency Correlation 
of Cavitation Surge with Rotating Cavitation and 
Synchronous - Water Flow 

Figure 1 1 overlays the normalized 
instantaneous frequency variation of a composite 
3(rc - N), developed from those of the rotating 
cavitation and synchronous components, with that 
of the actual cavitation surge component. Along 
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with the overlay are displayed the variations 
(normalized by mean frequency over the entire 
record) of the rc and N components used in the 
construction of the composite 3(rc - N). 
Generalized Hypercoherence processing was 
used to develop the micro-frequency variations 
shown in Figure 11.* The commercially-available 
signal processing package PC Signal, offered by 
Al Signal Research, provided the powerful 
diagnostic technique. The convincing match of the 
actual and reconstructed instantaneous 
frequencies of the cavitation surge oscillation 
indicates a high degree of synchronization 
between the rotating cavitation and the cavitation 
surge evident in the water flow circuit. 

Instantaneous frequency analysis was not 
necessary in making this same assessment in the 
LOX environment. Figure 12 is a waterfall plot of 
inducer inlet pressure taken from the last hotfire of 
the slightly deformed hardware shown in Figure 3. 
The frequency axis has been scaled to facilitate 
comparison of water flow data with that of the 
LOX. During the first portion of the test, the LOX 
pump experienced an excursion into rotating 
cavitation (start +52 through 55 seconds). The 
LOX pump was operating at 88% of rated flow 
during the transient and traversed the critical Nss 
-egion of the rotating - synchronous cavitation 
ooundary (Figure 6, -1 6 kNss). If the 
suDsyncnronous component were locKed into a 
frequency variation corresponding to 3(rc - N). it 
would be a mirror image of the (3N - rc) 
component about the rotating cavitation frequency 
trace. The subsynchronous component does not 
; ollow this variation during this excursion. 

However, this frequency correlation between the 
subsynchronous component and the synchronous 
and rotating cavitation disturbances increases 
toward the end of the transient. This varying 
degree of correlation, similar to that documented 
n the LE-7 experimental effort, is thought to 
-eDresent a reduced synchronization oetween tne 
ocanzed inducer disturbance (rotating cavitation) 
and the overall hydraulic system (via the 
suspected cavitation surge mode). Envelope 
analysis of both inducer inlet unsteady pressure 
and pump housing vibration from LOX testing also 
exhibited a well-defined periodicity at the 
subsyncnronous frequency. 


iiinua 

M. I >M '-•CaJCX*X m <• a«MO 

tLT imi * m*<a wu/m mm -<>+>- mu *• 

MV* • ■ «■/ t 



Figure 12. Water Plot of Inlet Fluctuating 
Pressure from Fastrac Turbopump Hotfire 

The spectral component labelled as f2 in 
Figure 4 showed little dependence on Nss and <t> in 
water flow test. The frequency decreased slightly 
f rom 4 Hz with increasing Nss. While the 3(rc - N) 
mode appeared to be most active in the vicinity of 
the inducer, the f2 mode was much more 
distributed across the water flow loop. 


III. Hydraulic System Stability Analysis 
Analytic Model Overview 

In an attempt to better understand observed 
system couplings” between the water flow LOX 
rnducer and test facility, a lumped parameter 
model of the MSFC Inducer Test Loop was 
developed. The analytic modei is purely hydraulic 
in that no structural participation is considered. 

The model is composed of simple mertance. 
resistance, and compliance elements as well as 
additional terms specific to the pump component. 
These terms include the mass storage term due to 
variation in o. or mass flow gam factor, and a 
oumD pressure gain factor. ~he pumo was 
modeled with two characteristic equations 
consistent with the method of Oppenheim and 
Rubin used in their Pogo stability analysis of liquid 
rocket engines. Even though the the selected 
pump representation could be viewed as limited in 
light of the four-terminal pump dynamic transfer 
function experimental efforts performed at the 
California Institute of Technology, the treatment 
was felt sufficient for a first-oraer analysis of water 
flow system stability. ’ 1 
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The completed model consisted of 
60 elements, which included lines, purely resistive 
losses (e.g. valves), accumulators, and the pump. 
The modeling of the hydraulic lines was consistent 
with that of Reference 7. Wave speed reduction 
due to fluid air content and pipe wall compliance 
were considered and incorporated into the line 
elements. Pump gain and resistance were 
developed from steady pump performance 
characterizations of the inducer-impeller in water, 
while pump inertance was estimated using an 
empirical equation from Fenwick. 9 

A strategy for setting pump compliance and 
mass flow gain factor in the analytic model was 
developed prior to anchoring the model to test 
data. First, dimensionless mass flow gain factor, 
M, was set by inlet cavitation number, a, and an 
assumed dimensionless pump compliance. K. 
using the recently published onset condition for 
both cavitation surge and rotating cavitation 
developed by Tsujimoto, Kamijo, and Brennen. 

M and K were converted to their dimensional 
forms. M' and C, respectively, in a manner 
consistent with scaling employed by California 
Institute of Technology researchers in their 
experimental determination of four-terminal oumD 
dynamic transfer functions. 

An initial estimate of K made using a maDPing 
or an estimated K and M versus inducer tip 
cavitation number, ay, published by Brennen and 
Acosta for the Rocketdyne J2 LOX turbopump. 
Even though the parameters were anchored to 
inlet perturbation testing where fluctuating mass 
flows were neglected, the normalized mass 
storage terms were deemed valuable because the 
J2 LOX turbopump was known to exhibit a strong 
supersynchronous oscillation thought to be 
rotating cavitation. 9 ,J The current shrouded 
version of the J2 LOX turbopump integrated with 
the Rocketdyne XRS-2000 AerosDike engine also 
exmbits rotating cavitation along with (rc - N). 


The experimental data attribute which drove 
analytic model adjustment of K (and subsequently 
M) was variation of the 3(rc - N) frequency 
observed with water flow suction pressure 

ramping. Figure 13 shows the reference 
frequency tracking taken from water flow test 
along with the last iteration of K(o.) and M (cr.). 
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Figure 13. K and M Estimates versus o T , Based 
on Observed Frequency Variation of 3(rc - N) 

System Stability Model Results 

A characteristic root of the system model with 
the pressure mode shape shown in Figure 14 was 
chosen to anchor K and M selection since it best 
matched observed amplitude distribution of the 
cavitation surge. Test data and audible 
observation showed the cavitation surge to be 
strongest in the vicinity of the inducer. The 
pressure mode, normalized by pump inlet 
pressure response, is plotted versus model 
element with flow direction going from left to right. 

Predicted critical damping for the moae ranged 
from -0.20 at 12 kNss to -0.52 at 26 kNss. 

Figure 1 3 also shows the frequency variation of 
the predicted mode versus a.. Interestingly, the 
most unstable modes predicted by the model 
consistently showed a mode shape similar to that 
of Figure 14, i.e. peak fluctuating pressure (and 
flow) at the inducer leading edge. 
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Figure 14. Predicted Mode Shape - Inducer 
Localized Instability 
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IV. Concluding Remarks 

Fastrac LOX inducer rotating cavitation shares 
many of the complex features observed in other 
rocket engine turbomachines which have exhibited 
the same phenomenon. 

The successful scaling of the oscillation between 
the water and LOX environments is encouraging. 

The unique interaction between rotating cavitation 
and cavitation surge at the inducer inlet warrants 
further investigation. A better understanding of the 
flow features which induce such complex dynamic 
signatures will most likely enable the design of 
very high suction performance rocket 
turbomachines in the future. 

The abundance of pump element assumptions 
made in the stability modeling attempt reflects the 
need for continuation of four-terminal, 
experimentally-derived, pump dynamic transfer 
functions as pioneered by California Institute of 
Technology researchers in the 1970’s. The 
severity of the cavitation-surge related damage to 
the Fastrac LOX inducer reinforces this need. 
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